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Abstract 

We show how coherent, spatially resolved spectroscopy can disentangle complex hybrid wave 
functions into wave functions of the individual emitters. This way, detailed information on the 
coupling of the individual emitters, not available in far-field spectroscopy, can be revealed. Here 
we propose a quantum state tomography protocol that relies on the ability to selectively excite each 
emitter individually by spatially localized pulses. Simulations of coupled semiconductor GaAs/InAs 
quantum dots using light fields available in current nanoplasmonics show, that undesired resonances 
can be removed from measured spectra. The method can be applied on a broad range of coupled 
emitters to study the internal coupling, including pigments in photosynthesis and artificial light 
harvesting. 

PACS numbers: 82.53.Mj,78.47.jh,78.67.Hc 
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I. INTRODUCTION 



The formation of collective optical resonances from Coulomb-coupled optical emitters 
is a very general phenomenon, including examples from chromophores in biological light 
harvesting complexes^', semiconductor quantum dots^, metal nanoparticles and composite 
systems, such as plasmon lasers^. 

For all these structures, dipole-dipole coupling occurs on a nanometer scale and the states 
of the individual emitters hybridize to form new collective, so called excitonic states, delo- 
calized over the whole structure. Far field excitation, governed by the wavelength resolution 
limit A/2, can only probe delocalized exciton states of a nanostructure. Related far-field 
experiments such as absorption, pump probe and four wave mixing are unable to disentan- 
gle the individual contributions of the coupled emitters from the collective optical response, 
because the exciting fields are spatially constant on the scale of the entire structure and 
cannot discriminate different emitters. In contrast, spatially local spectroscopy such as near 
field spectroscopy can, in principle, address the individual emitters. 

In this paper, we propose a new class of measurements that combine coherent nonlinear 
spectroscopy with near field optics to reconstruct the contributions of single emitters to the 
delocalized wave function in a spatially extended nanostructure. As an example, we demon- 
strate, how a coherent double-quantum-coherence optical technique^ may be combined with 
spatially localized fields to reconstruct the exciton wave functions of three dipole coupled 
self-organized GaAs/InAs quantum dots. This constitutes a particular quantum state to- 
mography. The presented procedure is independent of the technique for localizing the fields 
at individual emitters. Several localization methods are known and are already applied to a 
broad range of nanoemitting structures, e.g. using (metalized) near field fiber tips™ 3 , metal 
tip d 11 * 12 *, nano antennas^HHI an( j me tal structures combined with pulse shaped fields^^. 

Quantum state tomography is a development aimed at the direct reconstruction of wave 
functions or more generally the density matrix, first proposed by FancP^. The importance of 
quantum state tomography results from the fact, that the reconstruction and knowledge of 
the wave function opens the possibility to calculate new observables not related to optics at 
all. Examples include magnetic moments and transport properties. So far, wave functions 
are seldom directly accessible by experiments^. Recent advances include imaging of single 
orbitals using soft-x-ray pulses^ES an d the reconstruction of states^ES Applications so far 
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range from Spin 1/2 particles^, photon states using the Wigner functiorP^^, vibrational 
states^ to Josephson junctions^. In contrast to earlier approaches, the quantum state to- 
mography developed in this paper combines optical fields, highly localized in time and space 
with coherent 2D spectroscopy, using a sequence of light pulses with controlled envelopes 
and phases^HSSl 

II. EXCITONS IN COUPLED NANOSTRUCTURES 

As a typical example for coupled nanostructures with delocalized wave functions, we 
study three coupled self-organized semiconductor quantum dot d 6 1 33 1 34 1 , cp. Fig. [T^i). The 
quantum dot distance is assumed to be sufficiently large to have no electronic wave function 
overlap between the quantum dots. In this case we study interdot coupling in the form 
of dipole-dipole (or Forster) coupling known from selforganized GaAs/InAs quantum dots: 
Parameters like dot size, dot distances, coupling constants, and energy shifts are well known 
from theory^ and experiment!^. Each quantum dot is represented as a two level system. 
This is a valid assumption for quantum dots provided (i) quantum dots have no spin-orbit 
splitting and a big enough biexcitonic shift, (ii) are negatively charged or (iii) have spin-orbit 
coupling bigger than the inter quantum dot coupling d 37 * 38 l For selforganized quantum dots 
with sizes of 20nm and interdot distances around 40nm, the dipole coupling is about several 
/j,eV with a Lorentzian zero phonon line (ZPL) width of 7 = 1/ieV at low temperatures 
(e.g. T = 4A") 1 38 1 39 1 . We neglected the influence of the phonon side bands, since their 
amplitude in the spectra is one to two orders smaller than the amplitude of the zero phonon 
line resonance at low temperature d 38 * 39 4 

Three coupled quantum dots exhibit joint states: a ground state g, three single-exciton 
states ei, €2 and e% and three two-exciton states /1, /2 and ^3, cf. Fig. [T]d). The system has 
one triexciton state, but these states are of no relevance in a third order optical experiment, 
considered here. The ground state of the uncoupled quantum dots is not changed by the 
induced dipole-dipole coupling. The delocalized single-exciton states |e) resulting from the 
dipole-dipole interaction are composed of local, uncoupled quantum dot states \i) (quantum 
dot i in excited state): |e) = X}jC?|i). l e ) is an energy eigenstate of the coupled quantum 
dot system, cf the expansion coefficients. Similarly, two-exciton states |/) are composed of 
states with two local excitations at quantum dot i and j: |/) = J2i<j c ij\ij) ■ 111 general, 
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excited states of N coupled two level system emitters form a ground state g, N delocalized 
single-exciton states e and iV(iV — l)/2 delocalized two exciton states /. For our three 
dot case, we choose couplings between two quantum dots slightly stronger than to the 
third quantum dot (parameters given in Table [T]). Here, Hq includes along the diagonal the 
transition frequency local emitters modified by single and two exciton shifts, respectively. 
The offdiagonal elements describe interactions describing excitation transfer caused by e.g. 
dipole-dipole interactions. 

First, to characterize the system within far field spectroscopy, we calculate the linear 
absorption spectrum: 



k 12 



e 9 1 



(U ~ UJeg) 2 + 7 



(1) 



Here, \i eg is the dipole moment for ground state to single-exciton transition, u eg is the 
transition frequency and 7 the dephasing constant. 

The absorption spectrum of the coupled quantum dot structure is plotted in Fig. pf solid). 
The single-exciton states e±, overlap spectrally such that only e\ and ti are well 

resolved, contributes only with a spectral shoulder. Comparing coupled and uncou- 
pled(dashed) spectra, one recognizes, that the oscillator strength is originally evenly dis- 
tributed but strongly modified, since the dipole-dipole coupling forms excitons delocalized 
over the entire structure. 
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TABLE I. Hamiltonoperator in matrix form, a) The single exciton block and b) the two exciton 
states block. All values are given in fieV. The diagonal elements of the matrices are given as 
detuning to a mean gap frequency, a) uj gap for the single excitons, b) 2uj gap for the two excitons, 



with u. 



gap 



1.053eV and the two exciton shifts Vti = O.lfieV, Vt2 = —2.5fieV and Vt 



/:-! 



-1.5/xeV. 
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FIG. 1. a) Three dipole-dipole coupled self-organized InAs quantum dots, b) Exciton level scheme 
of the three coupled quantum dots. 
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FIG. 2. Absorption spectrum coupled (solid) and uncoupled (dashed) quantum dots. The detec- 
tion frequency ui is given as detuning relative to frequency Lu gap = 1.053eV (transition frequency 
of uncoupled quantum dot 3). 



III. INGREDIENTS FOR RECONSTRUCTING DELOCALIZED STATES 



Our main goal is to gain information on the built up of the delocalized wavefunctions of 
the excitonic states, i.e. on the expansion coefficients c\, for a given single- exciton state \e). 
For this purpose, we use coherent, spatially local spectroscopy, composed of three ingredients: 

(i) local nanoscale excitation provided by metallic nanoantennas and refined pulse shaping 
technique d 19 | 4C| to optically address individual quantum dots, (Section III A) 

(ii) phase cycling of the optical response^®', to disentangle the total nonlinear response 



into desired quantum paths, (Section III B ) 



(iii) a postprocessing procedure to calculate the coefficients c\ (Section IV). 



In general, (ii) and (iii) can be applied to any quantum system representable by spatial 
separated coupled emitters, if any localization technique (i) is available. 



5 



A. Localized excitation 



A main ingredient of our scheme is the local excitation of individual quantum dots. In our 
specific example, we achieve local excitation of the individual quantum dots by a plasmonic 
antenna structure of triangular symmetry on a subwavelength scale, cp. Fig. [3^,). These 
metal structures can be realized by e-beam lithography. Solving Maxwell's equations for 
this geometry shows that plasmonic effects and an optimization procedure of the applied 
pulses allows to selectively excite single quantum dot d 44 * 45 !: 

For optimizing the pulse envelope of a single pulse E(£, r) towards a field localization 
at only one quantum dot, we use time-harmonic solutions E„(u;, r), represented by incident 
plane waves of polarization directions p, s and incoming direction (indexed as 

1 f°° 

E(t, r) = / du £ 9u(oj)E v (u, r) e — *. (2) 
Pulse shaping is introduced by the weighting function: 

g v (u) = J2 /"( 1 ? n )^ e -("«-") 2 ^ 2 /2+ lW r„+^ ) (3) 
n V27T 

which represents a composition of Gaussian pulses with amplitudes A, center times r, fre- 
quencies T), widths a, phases f3, and polarization angle $ for each pulse n projected to 
polarization direction v (f p = cos, / s = sin). g u (u) has to be determined by optimization. 
To increase the number of optimization parameters, we combine the three incoming pulses 
from three directions, using 120° symmetry of the sample. For this paper, details of the 
optimization procedure are of no relevance but can be found in Ref. HQ] and HH Later on, 
the absolutes value of E(t) in the quantum dots centers is the input for the calculation of 
the localized spectra. 

In Fig. [3Jd) , the spatial field distribution for the optimized total field around the quantum 
dot transition frequency is shown. It can be recognized, that a chosen, single quantum dot 
is excited stronger than the other quantum dots. We observe field enhancements between 
different quantum dot sites of a factor of eight or larger. Note, that the optimized fields in 
frequency domain show that polarization and propagation phase effects cause localization 
and not a frequency based selection of different quantum dots. 

Note, that the presented localization scheme using excitation pads is just an example. For 
application of the protocol to other systems^HIS other spatial localization schemes might 
be used. 
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FIG. 3. a) Schematic geometry: Three selforganized GaAs/InAs quantum dots (diameter 20nm, in- 
ter dot distances 40nm) and three 12 nm thick silver layer structures, arranged with 120° rotational 
symmetry on a semi infinite GaAs- layer, b) Optimized localized electric field \E\ at the maximum 
peak for a single pulse composed from shaped pulses from three different directions. Note: the 
white color is 5 arbitary units or higher. The magnitude (extracted at X) of the optimized electric 
field \E(r ,tdgap)\ for ojgap at dots 2 and 3 is 12% or 7.5% compared to dot 1. The field intensity is 
scalable while third order perturbation theory is valid for the double quantum coherence spectrum. 
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FIG. 4. a) Pulse sequence for the double quantum coherence experiment, b) The two density 



matrix pathways Si and Su ( Eq. 20). See text for more details 



B. Phase cycling detection of coherent signals 



As explained in Sec. |III A\ a sequence of three spatially optimized pulse envelopes E l 
with phases (ft and laser frequency ui is usecP^S 



E(r,t) =E\v,t-b i -t 2 -t l )e 



iui(t-t 3 -t 2 -t 1 )+vp 1 
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FIG. 5. DQS for t% = 200ps a) absolute value b) imaginary part, J7i (O2) given as detuning 
around the single (double) gap frequency Lo gap . c) Imaginary part of localized double quantum 
coherence spectrum, where the first pulse is localized at quantum dot 1, 2 and 3 as indicated d) 
Filtered standard double quantum coherence spectrum (e2 removed). 

+E 2 (r, t-t 3 - t 2 ) e ^('-*3-*2)+^2 

+E 3 (r, t - t 3 ) e ^(*-'3)+^ + c.c. (4) 

Here the envelopes E l (r, t) are determined by the optimization procedure for localized pulses. 
The detected signal (selected quantum pathways of the full dipole density) is measured with 
heterodyne detection via phase cycling ^ 43 1 47 1 by repeating the experiment several times for 
different phases ipi, (f2 and ip$, cf. Fig. |4^i). 

In general the polarisation, created by three pulses applied to the quantum dots, is 
described by many quantum pathways in Liouville spaced In the following way, we can 
extract a subset of the Liouville pathways by extracting a certain phase combination of tpx, 
(P2 and (f3\ The detected dipole density for different phases can be written as^ 

P(t, (pi, (p 2 , ¥3) = V{t, cpt, (f 2 , <pa) + c.c, 

V(t, (fix, tp 2 , (£3) = C 123,lmnPlmn(t) , (5) 
Iran 

with c 12 3,imn = e z{l,pi+mip2+n ' P3) , / + m + n = 1 and \l\ + \m\ + \n\ = 1 or 3 for resonant 
excitation and Pimnif) being the part of the detected polarisation with phase dependence 
lifi + nup2 + n<f3- Ci23jmn can be viewed as a matrix with first index (ipi, (p%, ips) and 
second index (l,m,n). Carring out the experiment for sufficient phase combinations <px, 
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(f2, y?3, so that the matrix c\ 23 ^i mn is invertable, we can extract the signal with a specific 
phase combination <^ 4 = lepi + rmp 2 + n<p 3 (selecting particular pathways) using: Pi mn {t) = 
Si,2,3 c_1 i23,imn'P(^ <Pi, f2i fz)- Details of this phase cycling procedure can be found in Ref. 
HU Typical examples for such signals are the photon-echo v?4 = — ipi + ip 2 + <p 3 , anti-photon- 
echo (p4 = <f i — if2 + ^3 (cf. Ref. IHJ). 



C. Double quantum coherence signal 

We focus on the double quantum coherence signal, a third order signal with the contribut- 
ing phase combinations = i pi+ i p 2 — <pg^. In the case of a system, where the ground state, 
single exciton and two exciton states form three bands (cf. Fig. |l]b)), only the two Liouville 
pathways depicted in Fig. [4J3), will contribute to the signal with <^ 4 = (p% + <p 2 — v?aP^- 

In the case of the three band model (Fig. [lb)) only two Liouville pathways can contribute. 
The part of the polarization attributed to y? 4 , i.e. P[ i_i(£) which depends on the delay times 
can be written using a reponse functiorP: 

/•oo rco roo . . 

/ dr 3 / dr 2 / dr^W^t, t - r 3 , * - r 3 - r 2 , t - r 3 - r 2 - n). (6) 
Jo Jo Jo 

Note, that we include the optical fields into the definition of the response which is rather 

uncommon, but for the use of localized fields this notation will simplify the discussion. The 

(3) 

response function R[ l_ t can be divided into the contributions of two Liouville pathways, 
extracted from the full response functiorM 

R i3 \t, t 3 , t 2 , h) = (£) 3 tr(/i£(f - t 3 )H el ^(t 3 )g(t 3 - t 2 ) 

H el _ L ^{t 2 )g{t 2 - ?i)£r rf _i,_(fi)po). (7) 

Here, the electron electric field interaction Liouvillian H e i_i-{t)p = [H e i^L(t), p], the Green 
function Q{t) with Q{t)p{t) = 9(t)exp(—^Hot)p(t)exp(^H Q t) and the dipole operator p = 
Si gild) (A 1 + h.a.. For our excitonic three band system, for the far field excitation we insert 
the light matter Hamiltonian in local basis: 

H*- L = £/v E(t)\g)(i\ + Y.figi ■ E(t)\j)(ij\ + H.a.. (8) 

i ij 
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The Hamilton operator can also be reformulated in the delocalized basis: 

H e i- L = $>*e • E(t)\g)(e\ + $> e/ • E(t)\e)(f\ + H.a., (9) 

ef 

with the delocalized exciton dipole matrix elements fx ge = Y,i cff^gi an d A*e/ = 2i<j c !*AV c i(r 

I I I I /q\ 

We insert Eq. (9) into Eq. (7) and collect for R\(_i only the terms proportional to 
exp[z(— <£> 3 + <£>2 + y?i)] and end up with the response from two contributing Liouville pathways 
(Fig. Ep)) assuming no temporal pulse overlapP: 

R ( { i J_ 1 (t,fa,fa,fa) = Ri 3 (t,fa,fa,fa) + Ru (t, fa, fa, fa) (10) 
Rf\t,fa,fa,fa) = - jonfr+h-h-ti-**-*) 

^ He'fUge' 4 E 3 *(fa — fa)flf e • E 2 (fa — fa — fa)fJ> e g 4 E 1 (fa —fa—fa— fa) 
ee'f 

e -^/ e '(*-*3)-^/ 9 (t3-f2)-*fes(t2-fl) 

Rf(t,fa,fa,fa) = {^\m(h+h-h-t 1 -2h-t a ) 



^2 Vge'He'f 4 E 3 * (fa — fa)flf e ■ E 2 (fa — fa — t 2 )fi e g 4 E 1 (fa —fa—fa— fa) 

ee'f 

e -< £ i g {t--h)-i(,f g {h-t2)-i^ 3 {h-^) (12) 

Here, £ nm = u nm — i^ nm , with u nm = uj n — u m including the exciton frequencies u n and the 
dephasing/relaxation rate 7„ m for a Lorentzian dephasing model. 

In both pathways (i,ii), we have a coherence between the single exciton and ground state 
in between the first and second pulse and a two exciton to ground state coherence in between 
the second and third pulse. After the third pulse the system is either in a single exciton to 
two-exciton coherence (pathway (i)) or ground state to single-exciton coherence (pathway 
(ii)). We consider for further analysis the heterodyne detected signal, where the emitted 
signal Pj j is mixed with the field of a local oscillator E±: 

s { ffjfa,fa,fa) = r ttp^Um'^ty^ (is) 

J — oo 

Sj® (fa, fa, fa) is a complex quantity. A measurement obtains the real part of Sj^ (fa, fa, faf^- 
However the use of a local oscillator in heterodyne detection allows -by twisting its phase - to 
detect also the imaginary part of the signaP^ 31 | 48| E3 (phase cycled detection of fluorescence 
in fourth order^can give similar information as heterodyne detected signals in third order), 
this works both for the signal in temporal and Fourier domain. It is therefore a prefered 
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method to extract also the phase information of the coefficients cf , most other methods will 
only allow to extract the absolute value. 

In order to separate the different coherences of the signal by their energies, the signal is 
Fourier transformed over the delay tiniest 

5g I (n 1 ,n a) t 8 )= ( / o dfcjf dhe^+^S^Jt^h). (14) 
For the analysis, the double quantum coherence signal S^. is plotted as a function of the 



frequencies fl±, Q 2 , cp. Fig. m 

sff n (n u n 2 ,t 3 ) = s^in^M + ^^u^ts) (is) 



11 ee'f 



A*e/ " * (-/e)^e ' E (u eg ) ^ - ^ ^ _ ^ (16) 



3J 

= -p /V' ■ E 4 *(u e >g)fi e ,f ■ E 3 *(u fe/ ) 

n ee'f 

/v • s ("/ e )/v • * W(^-e /fl )(^-e e9 )- (17) 

It exhibits resonances for the ground state-single-exciton transitions u eg along the Qi 
axis and the ground state-two-exciton transition cu/J^S along the f2 2 axis. Due to the use 
of the local oscillator, the imaginary and real part of (Cli, Q 2 , £3) can be obtained from 
experimental dataP. 

D. Localized double quantum coherence signal 

For localized spectroscopy described here, the double quantum coherence signal S^p^ 
is modified by localizing the first pulse at a specific quantum dot i, cf. Fig. [3)3). 
For a localized excitation the Hamiltonian Eq. (pi) must be modified: 

H a - L ^Y,^9i- E {nM9){A+Y,^9i- E {nMj)Sj\i + H -^ (18) 



and yields 



H el _ L = Y,c^ gi -E(r l} t)\g)(e\+ ]T cf ^ ■ E(r i: t)\e)(f\ + H.a. (19) 

ie i<jef 
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for the delocalized states. We see that no delocalized dipole moments are formed, since the 
effective response depends on the spatial distribution of the electric field. 

Using far field excitation for pulses E 2 , E 3 , the local oscillator E 4 for heterodyne detection 
and a localized excitation for the first pulse E 1 at dot i (E 1 — > El), the double quantum 
coherence signal Sj® (i, Oi, 2 , t s ) = s!- 3 \i, 2 , t 3 ) + S^(i, Qi, 2 > £3) now dependents 
on the chosen quantum dot i and reads: 

sg^i, q u n 2 ,t 3 ) = sl 3 \i, n 1} n 2 ,h) + s£\i, n x , n 2 , t 3 ) (20) 
7(3)/ 



5r(i ) o 1 ,o 2 ,t 3 ) 
1 



3 S Me'/ ■ E *(^/e')/V ' -£ 3 *(^e' 9 ) 
ee'/i 



* rf/ \ e* * pi/ \ exp(-«f /e /t 3 ) 



(21) 



4 ?) (i,fii,fi2,t. 



3 

1 

h 



3 H /V • E 4 *(u e/g )fi e > f ■ E 3 *(u fe 



ee'fj 



Al • E* Mc> *„ • ^ (rjlW .)^^L. (22) 

fi e g/Hf e are single-exciton/two-exciton to ground state/single-exciton dipoles in the delocal- 
ized basis and \i g i is the dipole moment for the ground state to excited state transition of 
quantum dot i. E}{ij,u eg ) is the first pulse, predominantly exciting quantum dot i, only 
weakly exciting the other quantum dots with i ^ j. We assume ideal localization by taking 

E l(rj,UJeg) ~ SijElir^Ueg). 



E. Discussion of the double quantum coherence signal 



Fig. [5] shows the far-field double quantum coherence signal (El(rj,cu) ~ E 1 ^), cf. 
Sec. HI C[ ) absolute Fig. 4a) and imaginary value b): The frequency of the single exciton 



to ground state coherence can be seen on the Qi axis and of the two exciton to ground 
state coherence on the £l 2 axis. Clearly, for the far field excitation in Fig [5^) and b) we 
see resonances connecting to coherence of several states e and /. If we select a frequency 
Qi = u> ei g, we see along the Q 2 axis, which specific two exciton states are connected via 
dipole moments to the single exciton state e« and vice versa. A comparison of the dipole 
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moments connected to two different peaks works only roughly, since two Liouville paths 
interfere and the degree of destructive interference is different for every peak. 

A dominant peak (A) in the absolut value spectrum (Fig. |5p,)) is connected to 62 and 
/2, a second strong peak is connected to €2 and f\ and some further peaks with smaller 
oscillator strength can be seen at a lower single-exciton energy (e^ and f2,ei and fi). t\ and 
€2 are well resolved, e% shows up as a spectral shoulder. This shows that the system has 
three single-exciton and three two-exciton states. 

Fig. shows the signal with the first pulse localized at either quantum dot 1, 2 or 3. The 
localization of the first pulse gives information about the single exciton states contributing 
to the ground state-single exciton transition occuring during the first pulse. Localization 
at quantum dot 1 shows that all resonances connected to the delocalized exciton state 
e\ disappear. Overall, this shows, that quantum dot 1 only contributes strongly to the 
formation of single exciton state €2 and e^, but not to the build up of e\. Similar information 
is obtained for excitation of quantum dots 2 and 3 (see other Figs. (5}j)). E.g. the exciton 
state €2 is formed by quantum dot 1 and 2. Another interesting feature is the peak connecting 
e% and f\. This peak is only visible at the localized spectrum at QD 2 and QD 3 and not in 
the far field spectrum. This is caused by the fact, that is an antisymmetric delocalized 
state between QD 2 and 3, seen by the opposite sign of the peak in the QD 2 and QD 3 
spectrum. For far-field excitation, these two antiparallel dipole interfere destructively, so 
that the resonance is not observed. 

We next use the localized double quantum coherence to extract the wavefunction coeffi- 
cients c\ and therefore all quantum dot interactions. 



IV. EXTRACTING THE SINGLE EXCITON WAVEFUNCTION 



All ingredients are now available to extract the single exciton wavefunction. We start from 



the localized signal in Eq. (20 22 ) and see that the sum over e and j prevents us to extract 



a particular coefficient cf. Assuming ideal localization of the first pulse at a particular 
quantum dot i [E] {r j , uj eg ) ~ 5ijEl(ri,u eg )) removes the sum over j in Eq. (20 22). Of 



course, any deviation from ideal localization will result in an error in the measurement of 
the coefficients (see below). 

For removing the sum over e and selecting a particular single exciton state e, we choose 
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the frequencies Q± = VL\ and = VL\ in a way, that only a specific peak caused by single- 
exciton to ground state u eg and two-exciton to ground state coherences Uf g connected to e 



contributes, as suggested by the denominators in Eq. (20 22). Again, if peaks for different 
single exciton states overlap, errors are introduced to the reconstruction. (However two 
dimensional spectroscopy has less spectral overlap than one dimensional spectroscopy, since 
the peaks are separated by an additional degree of freedom: the additional frequency axis.) 
This yields: 

s^ji, fif, n e 2 ,t 3 ) = s?\i, ni, nit 3 ) + s$\i, n?, fi e 2 , t 3 ) (23) 
1 



e'f 

* T7i2/ \ e* * 771I/ \ exp(— l^fe'ts) , nA \ 
Hef ■ E (W/e)Cj fi gi ■ E i [n, U eg ) — — — (24) 

( i2 2 - Qfg)[ih - Keg) 

- -pE^' " E 4 *(u e/g )fi e/f ■ E 3 *(uj fe ,) 

11 e'f 

* 7-i2/ \ e* * itiI/ \ exp(— l^e'gt'i) / oe x 

Me/ • ^ (w/ e )Ci • (r f , w,*) ^ y — --. (25) 

l"2 — Kfg)\}h — Keg) 

We see, that here the double quantum coherence signal is proportional to cf*/i*j-£^ J 1 (r i , u; eg ),i.e. 
to the strength c e * the z-th quantum dot contributes to the delocalized wave function. This 
fact is used to develop a scheme to extract the coefficients cf from measured data: 

As input information the dipole moment \x g \ of the individual uncoupled quantum dots 
are required, the dipole moments can be measured or calculated. 

As measurement, carry out the localized double quantum coherence signal Sj® (i, fii, f2 2 , £3), 
for a localization on all quantum dots %. If the field strength and polarisation direction is 
different for localization at different quantum dots, we need to obtain the electric field along 
the local dipole \i* gi ■ E}(ri,cu eg ). 

Now, we select the excitonic state e = e a , whose coefficients c| a should be extracted. We 
determine the the frequencies fl± ~ u ea g, ^2 ~ Uf„ g showing a strong correlation to e a using 
the double quantum coherence signal without spatial localization. 

Now in the postprocessing of the data, we use that cf a * oc S^ n /{ii* gi ■ E}{vi,uj eg )) at 



the positions Vti « u ea g, ^2 ~ (Eq. (23 25)). c^ Q * can now be determined up to an 
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proportionality factor A: cf a *A = S kiii /(fi* gi ■ E}(ri, 0J eg )) for every quantum dot i, using the 
same frequencies Qi, fi 2 - Since the wavefunction is normalized, \A\ 2 = Ya l^ c i a *| 2 holds. We 
thus get A up to a global phase and set A = \A\. We obtain c e * = Sj^ n (i, Qi, f2 2 , t 3 )/(jj,* gi ■ 
El(Ti,u> e g))A). This gives the delocalized wavefunction \e a ) = J2i c i a \i)- 
Note, that these steps constitutes a quantum state tomography. The local basis is uniquely 
determined up to an arbitrary phase for every quantum dot: the expansion coefficient c\ 
depend on that choice. 



C.8 
0.6 
0.4 
0.2 


QD 1 QD 2 QD 3 

FIG. 6. Original and reconstructed coefficients of single-exciton wave function e 2 : Phase in 
multiples of 2ir. Error of absolute values determined by localization. 



To demonstrate the success of the tomography, we compare in Fig. [6] the elements of 
the reconstructed wavefunction for the strongest contribution, i.e. state e 2 (marked with 
A in Fig. [5^,)), to the original wave function resulting from the input parameters in the 
Hamiltonian. The agreement for both the amplitude and the relative phase is quite good. 
The difference results from a non-perfect localization E}{Tj,uj eg ) ^ 5ijE}(ri,u eg ) resulting 



from realistic Maxwell simulation from section |III A[ This error is marked by the error bars 
in Fig. [6j It is caused by a weak excitation of quantum dots, which a ideally localized pulse 
should not excite. Such a non ideal excitation leads to a cross talk between the coefficients. 
The error bars are estimated to be smaller than: A c e = \Ei{ T j)\/\Ei( T i)\- 

Note, that in general, the procedure works also for other methods than heterodyne de- 
tection in far field, including a localized detection of polarisation or fluorescence, as long as 
the detection is the same for a localization of the first pulse at different quantum dots. The 
only limitation is, that the phase of the coefficients can only be detected with methods, that 
can measure complex signals. For other types of detection like homodyne detection, we can 
also extract the absolute value of the coefficients, but not their phase. 
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V. FILTERING COHERENT SPECTRA 



As additional useful application, we show that strong, undesired resonances can be selec- 
tively suppressed from coherent spectra. This can be advantageous while investigating weak 
resonances, that are masked by other strong resonances: Often, it is not clear, whether weak 
resonances constitute a vibrational side peak connected to a dominanting strong excitonic 
peak or a different, much weaker excitonic resonance. This can also be solved by selectively 
removing excitonic resonances from measured spectra, applying a filter algorithm. 
As input information for the filter algorithm, we have to determine expansion coefficients 
c\ a for the specific state e a for all quantum dots i, whose contributions we want to filter out. 
Additionally, we need all dipole moments \i gi of the individual nanostructure and also the 
electric field along the local dipole \i* gi ■ E}(ri,u eg ). 

As measurement we record a localized version of the spectrum to be filtered. The localized 
pulse should excite a ground state to single exciton transition for all quantum dot positions. 
For the double quantum coherence, this will be the signal Sk IZI (i, ^1,^2, £3) for every quan- 
tum dot i. 

For postprocessing we discuss the expression 

Sk IU ,w/oe a (^1, ^2) £3) 

i 

= £ 4*^ g -El(u eg )fA^i^2,t 3 ) (26) 
f ea (Q l7 Q 2 , t 3 ) = £ d?S kl „ (i, Oi, fla, t 3 )/(/4, • El(r t , u eg )) (27) 

i 

which gives a spectrum, where all contribution of e a during the first pulse are filtered out.^ 
The single-exciton peak e 2 (a = 2) dominating the spectrum in Fig. [5] a) is filtered out 
in Fig. pp.) . This spectrum reveals now information about states initially covered by the 
dominant contribution of e 2 . The procedure can be applied iteratively, using the filtered 
spectra for obtaining the other excitonic states. This can enhance the reconstruction of the 
exciton states. 

The filtering method can also be applied to other spectroscopic signals as long as a phase 
sensible detection is used and a localized signal, whose contributions are proportional to the 
single exciton expansion coefficients, can be measured. 
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VI. CONCLUSION AND OUTLOOK 



The presented quantum state tomography for the extraction of the delocalized single 
exciton wave function coefficients, can also be applied to other impulsive two dimensional 
spectra. The single-exciton to two-exciton transition in double quantum coherence using the 
localization of the second pulse also also can be used to extract the two exciton coefficients. 
However since this problem is more complex, it will be subject to future work. 

In conclusion, our simulations demonstrate a quantum state tomography that can be 
used to reconstruct individual wave functions of coupled emitters acting only collectively 
in the far field. In addition, localized excitations are useful to remove unwanted strong 
resonances to uncover weak or hidden excitonic resonances. All of these features are not 
accessible in standard far field spectroscopy. Similar configurations can be alternatively 
achieved by applying four pulses and using phase cycling to detect a desired component 
e.g. with phase ip = <fi + (p2 — y?3 — fi- We therefore believe that the proposed quantum state 
tomography opens a new path for the detection of many body interactions on the nanoscale. 
The proposed protocol is more general as presented here, since fluorescence can also be used 
rather than heterodyne detection of optical fields®!. 
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s klII (i,ni,n 2 ,t 3 ) = Y,4*i4g-EhTi,ue g )fe>(n 1 ,n 2 ,t 3 ) (28) 

e' 

Multipling the equation with cf and summing over i yields a scalar product and we get f e i 
defined using the localized spectra (Eq. p7[)). 

The far field double quantum coherence spectrum can also be calculated using f e r. 

s klII (nx,shM = Y,^'* ^ 9 ■ El Mfe'(^i,n 2 ,t 3 ). (29) 

e' 

It is expressed using summands for every contributing single exciton e', the summand of the 
exciton to be filtered can be substracted. This is possible, since every summand can be calculated 
using / e ', which can be calculated from the experimental data of the localized spectrum, if the 
expansion coefficients of the single exciton wavefunction of e' are known. 
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